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Abstract 

Using for unperturbed electron and phonon Hamiltonians a representation by the Jacobi 
matrices a one-dimensional model of the electron-phonon interaction is constructed. In frame 
of the model the polaron and scattering spectral bands are calculated explicitly as well as 
the corresponding states. An explicit formula for the polaron effective mass is derived and 
the fit of model parameters is suggested. A connection of the model with the description of 
phonon degrees of freedom for quantum wires is discussed. 
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1 Introduction 



The physical motivation of the present paper is that even for quantum description of the 
simplest nanoelectronic devices (like a quantum wire) one have to take into account the 
interaction of the propagating electron with the internal degrees of freedom of the device. 
This interaction becomes valuable in the case for instance, when the voltage applied is large 
enough or the temperature of the device is high. The simplest excitations in both cases 
are phonons. It means that the first step in the correct description of any nanoelectronic 
device at nonzero temperature or at large voltage applied is taking into account the electron 
phonon-interaction. This leads to polarons as charge carriers instead of electrons. In general 
the effective mass of polaron is greater than the one of electron in the lattice. So one of the 
important problem is to calculate the effective mass of charge carriers and to estimate its 
contribution into the physical characteristics of the working regime of the device. 

The standard field theoretical models the electron-phonon interaction such as Frolich 
model [1] do not preserve the number of particles and cannot be solved explicitly. This makes 
interesting the construction of the phenomenological models of electron-phonon interaction 
with fixed number of particles. The first step in this direction is the study of two-particle 
models (one electron interacts with one phonon). 

The first attempt to construct such a model was undertaken in the paper [2] where 
the electron was described by the one-dimensional Schrodinger operator with the Kronig- 
Penny potential and phonon was described by the Jacobi matrix. The interaction in this 
model is localized on the diagonal of electron-phonon space. Being exactly solvable this 
model however is not able to calculate explicitly the effective masses of the quasiparticle 
aroused. The physical difficulty related to this model is the difference of scales introduced 
for description of phonon and electron (continuous configuration space for electron and the 
lattice for phonon). 

In the present paper we describe and solve a similar model where the kinetic energy 
operators are given by Jacobi matrices for both electron and phonon. This makes possible to 
calculate explicitly the eigenfunctions of the polaron branch of spectrum and corresponding 
effective masses. Another simplification in comparison to the model from [2] is that we take 
into account the influence of the lattice potential on the electron by changing of its mass by 
the effective one in the kinetic energy operator. 

2 Model 

We start with the one-dimensional lattice Z. Set H = l 2 (Z) be the state space of a quantum 
particle on the lattice and T be the unitary shift operator in H, (Tip) n = tp n -i. Introduce 
the unperturbed phonon and electron hamiltonians on H as L\ = T + T* and L 2 = £(T + T*) 
respectively. Here the real parameter £ makes sense of ratio £ ~ j^£o_ U p ^ a di meri sionless 
factor where u is the frequency of the lattice oscillations and E c h ar is a characteristic scale of 
the electron energy i.e. the width of conductivity zone. The hamiltonian of electron-phonon 
system without interaction acts in 7i = H ® H = / 2 (Z 2 ) as L = L 1 <g> / + / <g> L 2 . Its 
spectrum is given by a(L ) = a{L x ) + a{L 2 ) = [-2(1 + 2(1 + 

We will suppose that the lattice is homogeneous and the interaction between electron 
and phonon is localized. This restricts the choice of the interaction operator V, L = L + V 
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being the perturbed hamiltonian by the following assumptions. 
Assumption 1. . 

[V, R] = 0, where {Rg) m ,n — Sm+i.n+i is the operator of the electron-phonon shift in H, - 
space homogeneity. 
Assumption 2. 

VnG Z g nn = implies Vg = 0, - locality. 

The locality assumption means that V differs from zero on the main diagonal V = {g e 
Z 2 (Z 2 ) : g mn = for m ^ n} of / 2 (Z 2 ) only. By the assumption 1 it commutes on V with 
the shift operator. Thus V is given just by convolution with a function \I> e Z°°(Z), 

n,rrt 0~n,m ^ ' ^ r9n~T,n—T' 

rez 

Assuming the symmetricity of V one gets = \I/_ r . If \I/ T is well-behaved (for instance 
if J2 T \ ^t\ < oo and so V is bounded) then V is selfadjoint. Generally one can define the 
selfadjoint operator V in T> ~ ^ 2 (Z) as the Fourier transform of the operator of multiplication 
on real valued function = J2 T ^t^ pt in -^ 2 [ _ ^, tt]- 

2.1 Spectrum of the Operator L 

To study operator L it is convenient to use the spectral representation of R given by the 
Gelfand transform T : / 2 (Z 2 ) -> L 2 ([-ff, tt], / 2 (Z)) = W, 

(^).(«) = 4=E^,ne* W . 

Properties of T 

1°. is isometric surjection. 

Proof 

Set D s = {g n>m 7^ =>- n — m = s} hence D s _L D s / for s ^ s' and 7i = ® st =zD s . Denote 
as D s = JFD S = L 2 (— 7r, 7r)e s , e s being the s-th element of the standard basis in / 2 (Z). We 
see that D s _L £> s / for s/s' and H = ®D S . Then JF| £)s is just usual Fourier transform. 
Thus J-'Id.s is surjective isometry. • 

2°. T diagonahzes R, R = F~ l e iK ■ T. 

3°. L = T~ X L§T where L = @L {k)oIk and 

L (K) = ^ + e'^)T+(C + e tK )T* 

We remind that T is the shift operator in / 2 (Z). 
4°. T diagonahzes V, V = T~ l \ ' T ', where 

(V^),(«) = <S i)0 tf («)&(«)• 

Proof 

Since -D s are invariant subspaces of R and JF the property 2° expresses the fact that 
shift turns to multiplication on function e lK under the Fourier transform. This together with 
commutation of T with the shift on s-variable implies 3°. The property 4° expresses the 
fact that under the Fourier transform the convolution becomes operator of multiplication on 
function. • 
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Thus L = Lq + V is partially diagonalized by T i.e. 

L— @L(n)dK, 

J —IT 

L(k) = (f + e- lK )T + (£ + e* K )T* + ^( K )P 

where (Po^Os = ^s.o^s is the projection in the fibre / 2 (Z). 

Since is the rank one perturbation of the standard Jacobi matrix we have a c (L K ) = 
[—2 |£ + e m | , 2 |£ + e m |]. This spectrum corresponds to scattered waves type eigenf unctions 
of L. It fills the zone 

<Jsc.wave= IJ °c (1 («)) = ["2 (1 + , 2 (1 + 

which coincides with the spectrum of unperturbed operator L . 



2.2 Polar on 

We are now interested in the waveguide type eigenfunctions of continuous spectrum of L, i.e. 
those which satisfy the condition <p n ,m = O (exp(— C\n — m\)). In physical literature such 
eigenfunctions are called polarons. It follows from the previous section that polarons with 
definite lengthwise momentum k G [— 7r, it] i.e. satisfying the condition Rip = x<P, X — e * K 
correspond to / 2 -eigenfunctions of L(k). This reduces the spectral problem L(k)<p = \<p to 
the one for Jacobi matrix 



/ 



\ 



V 



o c + x o 

o t+x o c+x o 

o c+x o e+x 

o c + x o 
o 



or 

Note that 

where 

and 



{Q x lf) n = (f + x)Vn-l + (f + XVn+1, ™ ^ 
(Q x ^)o = (e + X)<f-1 + + (£ + X)<Pl- 

U* k L(k)U k = \Z + X \ (j + v) 
U K = diag {e mT } neZ , r = arg(£ + x), Aj = s *j+i + s *j-i 



IC + xl 



Thus the / 2 -eigenfunction of Q x corresponds to the eigenvalue A=|£ + x|(a + a x ) given by 
solution of the secular equation 



■a - 1 = 
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with \a\ < 1. Solving this equation one gets that the polaron energy A is related to k through 
the dispersion relation 

A(k) = sign^(^) v /^(/t) 2 + 4(l + £ 2 + 2£cos K ). ( 2 ) 

The corresponding eigenfunction of L(k) has the form 

ip n (K,X) = e inT aW. 

Note that equation (2) shows that in the case of strong coupling 3> 1) the polaron zone 

is separated from spectrum corresponding to the scattered waves. For the small coupling 
^(k) ~ the polaron zone and a sc , wave are overlapped. Note however that the polaron decay 
is impossible for the energies belonging to the common part of spectrum. This is the result 
of high common symmetry of L and L$ provided by R. 

For completeness we now right down the generalized eigenf unctions of L(k) scattered 
waves type. The solution of the equation L(n)(p = \<p corresponding to A = 2|£ + x\ cos # is 
given by 

<Pn = • ' 

where 



S K (X) = 



e *M + S K {\)e' l9 \ n \ 
A-tt(/c)-2|£ + x|e 



A-*(«)-2|e + x|e- ifl 

Here 9 e [— n, it] is the quasimomentum corresponding to direction transverse to the diagonal. 
The reflection coefficient S K has the pole on the physical sheet in A- variable which corresponds 
to eigenvalue (2). Since Q x is the rank 1 perturbation of standard Jacobi matrix this pole 
is the only singularity of S K on the Riemann surface of the energy A. 



2.3 Effective Masses and Fit of the Parameters 



Using (2) one can now calculate explicitly the effective masses at the ends of the polaron 



zone by the formula m = 



\"(k) 



K = KQ 



provided A'(k ) = 0. The result is 



^(k) 2 + 4(1 + ^ 2 + 2£ cos k) 
\fr> 2 + ijn]/" _ 4£ C os k 



K=KQ 



where k is defined by the equation 

*(/Co)*'(ko) =4^sinK - 

The case ^>(k) = const or \1/ T = ip5 Tj o comprises a precise discrete analog of the interaction 
suggested in the paper [2]. The corresponding polaron zone is a, 



pol 



vV 2 + (iei-i) 2 , 



' 2 + (|C| + l) 2 for ip > and — o voi for ip < 0. The coupling constant ip here is an analog 
of the Frolich constant a [1]. This constant is incorporated in the Frolich hamiltonian H as 
H = H + ^faV ', where H is the operator of kinetic energy of electron-phonon system and 
V is the operator of electron-phonon interaction. In the limit of weak coupling ip ~ by 
(2) we have A — Aq ~ 4 | x ^p V ;2 ) where Aq = 2|% + ^| is the rest of the zone corresponding 
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to scattered waves. Comparing this with the asymptotic AA ~ a of the polaron energy 
measured from the rest of the continuous spectrum known from physical literature [3] one 
takes ip = Cy/a, C = 2\x + £|. 

Then we have sin k = what gives k = it for minimum. Thus 

vV 2 + 4(i-0 2 

m = — . (4) 

At ijj ~ the polaron effective mass (4) should coincide with the electron effective mass 
which is 1 in the accepted system of units. This condition fixes the parameter £ as £ = |. 
This defines completely C at the rests of the polaron zone. 

Let us note that the Frolich constant a can be expressed in terms of the physical char- 
acteristics of the material [4] 



a 



1 / , ,x e 2 ( 2m*u y 



where e and are dielectric constants at zero and high frequencies, m* is the electron 
effective mass. Thus we conclude that the coupling constant tp is completely fixed by the 
choice of the material. As both of the parameters ip and £ are fitted one can calculate the 
polaron effective mass (4) for various materials. 

In the weak coupling limit ip ~ electrons practically do not interact with the lattice 
phonons. Vice versa, in the strong coupling limit ip — > oo electrons strongly interact with the 
lattice phonons. In this case the description of the electron-phonon interaction in the frame 
of field theory becomes essentially nonlinear. This means that the linear approximation we 
used is not suitable in this regime. That is the reason why the asymptotics A ~ y/a, m ~ y/a 
at ip — > oo obtained from (2) under the accepted fit of ip differ from the physically expected 
A ~ a 2 , m ~ a A [5]. 

Note that when £ = the polaron zone degenerates into the point A = y/if) 2 + 4. This 
point is the infinitely degenerated bound state of L. One can interpret this as the result of 
braking the electron by vacuum of infinitely heavy phonons. 



3 Conclusion 

The suggested model of electron-phonon interaction demonstrates on qualitative level the 
rising of polaron branch in the spectrum of the Hamiltonian under consideration while the 
spectrum of scattered waves remains the same as for unperturbed problem. This model 
leads to the explicit formula for the polaron effective mass. The latter makes possible to fix 
the model parameters ip and £ as it is described above. Being applied to the description of 
the electron propagation through a quantum wire in non-ballistic regime the model allows 
to calculate the effective mass of charge carriers in terms of tfj and £. It opens the way for 
estimation of the contribution of electron-phonon interaction into the conductivity in the 
working regime of any nanoelectronic device based on quantum wires. The simplest such 
device is the quantum wire with the periodic external voltage applied [6, 7, 8, 12, 13, 11, 12, 
13, 14, 15, 16, 17]. As it is has been shown [6, 7, 8, 12, 13] in the effective mass of charge 
carriers is important characteristic of devices based on Mott-Peierls stimulated transition 
and could change the upper limits for the temperature and gate voltage applied. 
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